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BSLO&  fHB  SUB  FACE  OF  A  COMPRESSIBLE  FLUID 
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,Davlenlya  V  Glut*  Ssrtxinaemoy  .Zhid¬ 
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Assume  that  pressure  is  exerted  at  a  certain  point  I 

i 


0  on  the  free  surface  of  a  compressible  fluid.  This  is  ! 
then  propagated  along  this  surface  in  the  form  of  a  shock  j 
wave  (Fig.  l) .  We  will  assume  that  R(t)  describes  the  j 
propagation  of  the  wave  front  and  the  distribution  of .  ■ 
pressure' on  the  surface  and  is  a  known  but  arbitrary  ! 


]  function  of  time.  The  pressure  distribution  in  the  depths! 
|  of  the  fluid  Is  then  found  with  the  aid  of  a  quadrature  * 


f 

f 


Fig.  1 

j  1.  Pressure  distribution,  along  a  head  wave. 


h  cases  where  the  velocity  of  the  front  is  greater| 


I  than  the  velocity  of  sound  in  the  fluid,  the  disturbed 
!  region  will  be  bounded  by  the  characteristic  arc  (the 


!  head  wave)  AB  and.  a  circle  of  radius  at  (Pig,  l)  *  Let  us  f 

t  f 

I  ! 

|  select  an  axis  Or*  In  the  plane  of  the  free  surface  and  an  \ 

I  axis  Oz  perpendicular  to  Or. 


In  the  region  ACB  the  pressure  distribution  takes  j 


the  form} 


p(rv  t,  <)=-  U - 

0  K  J  J  (s*  +  r*  +  /$ — 2rf0  cos 

n 

jv  f*5  :  Pi  (ri>  *#) 

I  ^  ft  a  j  xs  +  r|  -t-  rp^j^fcosty  ^  ~~ 


<h  rg  gfefcJl  gf  fo  M 

,  j?  {  b%  y  » 

'  «ca  \  j  Xs  4-  f *  +  i*| — 2fr0  cob  4' 


*5  j  s8  + /f  +  rf*  cos* 


d,l) 


where  r^  g  are  the  roots  of  the  equation 


r*  %  =  R  ( t  —  8 1  j.  c0^  £ 


(1,2) 


|  f(r}  is  the  inverse  function  of  R(t),  Is  the  value 

|  of  ’f  for  which  r*(  ^)  =  r*(  ^  =  J*.  Then  and 

I  * 

\  f^  are  found  from  the  equations 


r*  „  ft  ^  __  Vj 'o* ± yg  ~  gVg, CQ3 ft  +  i!  ^ 


+i  J?'(L- 
■  «  v  ■ 


V /?*  +•  r|  ~~  2>jr0  cos  4>i  +  2s 


\ _ yg— y0cos  _ 

/  1/3$  +  bg*  ~  Syg  ‘ cos  +1* 


0.  (4,3) 


$  i 

|  The  head  wave  may  be  expressed  par ame t r i c a 1 ly  as  the  | 

I  envelope  of  the  elementary  disturbances  arising  at  a  I 


3 


point  on  the  free  surface  at  the  instant  of  time  t1  at 
which  the  wave  front  reaches  it, 

R'  (£')  [/?  (£')  —  ae  {tf  —*£).  (1,4) 


I  The  Eqs.  (1.3)  and  (1.4)  constitute  a  system  of  equations 

I  ^  *  *  • 

|  for  the  determination  of  rQ,  *px,  t'»  and.  .z  as  functions 
\  of  xv  and  t.  It  is  obvious  that  in  this  case  it  is  poss- 

i  0 


I  ible  to  set  ^ x  ~  0  and 


r  *< 


Y  »'e  +  ~  *'o  +  i1 


|  Then  the  Eos.  (1.4)  will  follow  from  (1.3)  and  for  the 

I  .  * 

I  determination  of  rri  and  z  there  will  remain  the  'two 

1  ’ 

Eqs.  (i«3) *  •  Thus  we  have  ascertained  that  ^  =  0  on 

"  AB.  Substituting'  U1  *  0  In  Eqs.  (1.3)  we  obtain  equa¬ 
tions  for  Tq  and. the  equation  of  the  line  AB.  Let  us 
j  denote  Tq  '  on  the -head  wave  as  Fq*  .It  la  clear,  that  this 
|  quantity  may  be  found  from  the  equations 


I  which 


L _ _  *  j 

Y(rS  —  /•«)*  +  «*  1  r*  —  r0 


**  I  —  ji 

J  VO* 


(r* — ,r0)*  f 


(1.5) 


which  serve  at  the  same  time  for  the  determination,  of  the 
equation  for  AB:  a  a  s(rr>s  t).  In  particular,  if  R*  s 
const  -  Vq  j  then  ■  "  . 
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t  *•  ^.«*-**-  ifc*  'fr'-^A'T 

t;  . •  ^ «-*>. - ^**4  •««•  /.*  iv  * »&.<**«&  ?  *•* *  ■HL*jt ?&<&.' fr n /<r-.‘jij': ft h W ■.\>’?ZWW^ ' Tfr-T57^; tf *s * ,: ■%'”■  ''*^<  a? ';^;. *  v  ?.  E  w •  ***•*•-  ^  ^ 

"d^e  iH'isg  Trt/ft.  After  substituting  this  value  of  r«  in 

t  or  Its.  (1.5) ,  w©  get  the  equation  of  the  head 

Ip^lfFroat  $ 


z  ~-----.ZZP.<!  * 
\At*~  i 


•:  Let  us  now  find  from  Eq,  (1.2)  an  expression  for  diq  g/dt, * 

l  ?  i 

|  .which  appears  under  the  integral  sign  in  Formula  ( 1*1) s  \ 


dr 


1,  $ 


dl 


14  --R'(t 
a  V 


w  ^  t  _ — ^  — 3^  **  ^  r  je*\ 

/  "£F 

*  *“  I.  VM-Jv  M  *t  -f  -  i  “O.  'll. 

i*  a 


V  *"*.‘a  4  r%~~  2Vj*  $  cos  41  -I-  a* 


r9  cos  4< 


V  r*\  4  4 — 2r<fl  J  cos  ty  +  z*  ' 

(m 


|  From  Eqs.  (1.6)  and  (1.3)  it  is  clear  that  at  the  upper 
j  ,  ■¥■ 

|  limit  (  4'  “  4*  \ )  ^rx  />/dt  becomes  infinite, 
i  ’  All  the  integrals  in  (l.l) ,  with  the  expect Ion  of 

i  * 

5 

j  the  last  two.  vanish  on  the  head  wave  AB.  In  the  two 

j  exceptions  the  interval  of  the  integration  degenerates 

}  to  a  point  while  the  integrand  goes  to  infinity.  From 

|  this  it  follows  from  among  other  things  that  in  the  case 

|  where  the  boundary  condition  is  given  in  the  form: 

1  Pi  (r>  0  «  PA  (#)  fi  C r}R ) , 

\  where  R(t)  is  the  radius  of  the  front  at  the  time  t, 

{ 

!  then  the  following  equations  will  hold: 

i 

1 


Pitt*  Q^PaK), 
Pi  (ri>  o  ~~  p&  {?*) 


t ■sm*'S*ru»  •■•* 

I 


that  la,  the  pressure  distribution  p(rQ,  z,  t)  along 
the  head  wave  AB  does  not  depend  upon  the  pressure  dls 


{  tributlon  behind  the  front  f-^r/R). 


;  -r  *0 


ffitt  Mr  SKW  mu  m  #»*»  *»»* «»<»  **  «**  «* 


8 

It 

» 

t 

i 

5 

4 

I 

I 

* 

I 


Fift-,  2 


I  To  determine  the  pressure  distribution  over  the  head  j 

5  I 

|  wave.  In  the  general  case,  we  must  find  the  limit  ap~  | 

|  proached  by  the  last  two  integrals  in  Formula  (l«l)  as  | 

i  we  move  to  AB.  ¥e  note  that  the  entire  continuous  part  J 

1  I 

of .the  integrand  may  be  taken  from  under  the  Integral  j 

sign  after  setting  •"  ? 

j  §**0,  7#.  j 

I.  i 

I  In  it.  .  Thus  the  problem  ultimately  becomes  one  of  finding 

\  '  f 

I  the  limit  of  the  expression i  s 


J  _Pi  (r$,  h)  ~  f  d  .  .  ,v  ,, 

««  *•  r0)*  0  \dt  (r* 


(i»7) 
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where'  7q  and  z  -  z  (rQj,t)  are  determined  by  Eqs.  (1.5)  •  j 
It  is  clear  from  (1*6)  and  (1.3)  that  only  the  function  ! 

a 

R(t*)  itself  and  its  first  derivative  HJ ( t; 1 )  enter  into  j 
the  Integral  in . Formula  (1*7)*  Aside  from  this,  the  value 3 

jj 

of  the  arguments  ;of  the  functions  R  and  R*  in  (1*7)  &P-* '  j 

*  i 

proach  t*  as  0,  that  is*  for  every  rQ*  2  and  t  I 

belonging  to  AB.  ! 

Let  us  consider  the  graph  of  R(t*)  (Pig.  2).  From! 

Eq»  (1.5)  we  find  and  t*.»  for  fixed  values  of  rQ,  a,  j 

s 

1 

and  t  on  AB*  •  | 

In  the  neighborhood  of  r  =  rA  and  t*  «  T*  we  sub-  J 

y  I 

stitute  the  equation  of  the  tangent  at  the  point  M  for  j 

% 

the  curve  r  -  R( t  • )  •  This  equation  takes  the  fora;  j 

r-rt-R'CW-Q  j 

The  eou&tion  for  the  determination  of  rv  0  now  1 


!  takes  the  form; 


>>*  ®s  4-  r| — 2 r®r*  t  cos  tf>  +  z* 


"  )  • 


¥e  obtain  an  expression  for  rj"#2  by  solving  this  equation: 


Af*!’®  cos  — f j  :fc  / lywfVj  cos  ’I' — fj)* — 

M* — 1 


(r|  4-  «9) 


(l.S) 


|  where 


»  r* «»  r?  -fs  M  (£#)  {t 1  i%)> 


We  have  the  following  expression  for  ^  : 

n  -f  |/  (M*  -i)  [Af*  (r|  +  ^sjl 


arc  cos  - 


iJ 


*1  —  ,  **« 

At  V0 

i  When  differentiating  (1.8)  over  time  it  is  necessary  to 
|  take  into  account  the  dependence  of  T*  on  t.  In  per¬ 
forming  this  differentation  it.  is  not  difficult  to  see 
|  that  the  integrand  will  consist  of  some  continuous 

!  i 

;  function,  which  we  will  discard,  plus  a,  fraction.  The  I 

i 

% 

numerator  of  this  fraction  is  also  a  continuous  function,  * 

\  1 

A 

I  hut  the  denominator  is  a  radical  which  also  enters  into 

f 

;  Eq,  (1.8),  and  which  vanishes  when  'f  =  ^  .  As  has 

i  .  ■  ■  I 

|  been  previously  noted,  the  continuous  part  of  the  inte-  I 
)  grand  may  be  taken  from  under  the  integral  sign  on  setting 
0  in  it.  Thus  the  entire  matter  reduces  to  the  cal-  ’ 


a,|' 


culatlon  of  the  integrals 
% 


V  cos  —r,j*  --  [M*  (rf  +  **>~  rfj  (M*— 1)  * 

i 

|  Substituting  in  the  last  integral  cos  tL'  as 

!  is  possible  because  of  the  smallness  of  ,  and  cal- 
!  culatlng  the  resulting  tabular  integral,  we  finally 
|  have  along  ABj 


AM 
dl 

7J®~~  i 


where 


n  (r  r  A  —  y~Mtra — n  I  rfrj  M  ~di  ^ 

F  {r0,  z,  1}  -  - -*»-*— 


(M> 
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and  the  pressure  distribution  along  the  head  wave  is 


given  in  terms  of  the  elementary  functions, 

i 

0,  rA  ~  R(t),  and  p(rn,z,t,fl 


On  the  free  surface  z 


~  p^( R, t) ,  If  Pj  =  const  and  R’ (t)  ~  const.  -  ?Q,  we 
have  for  the  pressure  distribution  along  AB t 


P  (ro*  0  = 


f M%  ~~V&t 


rs(M*~  i)  * 


(140)1 


The  Expression  (1.10)  was  derived  by  A.  Ya»  Sagomonyan  . 
by  the  method  of  characteristics® 

2,  An  Investigation  of  the . Solution 
in  a  Half  Space 

It  is  interesting  to  investigate  various  parti¬ 
cular  cases  of  the  definition  of  the  law  giving  the 

pressure  distribution  on  the  surface  of  the  fluid.  For 

* 

simplicity  we  will  assume  that  the  pressure  distribution 
to  be  uniform  on  the  liquid  surface  behind  the  front, 
but  to  vary  in  time;  p^(r,t)  =  p^( R) ,  where  R  «  R( t)  is 
the  equation  of  motion  of  the  front  over  the  surface. 

The  simplest  form  of  pressare  distribution  along 

the,  axis  0z  below  the  surface  will  be 

«* . . . . 

*ts 


up.  io- 


3  y  +■ 

© 


where  . 


,  VH*' +•*  ^ 

i  — j  j 


v?~+  ** 


v  ^  ■  i 

l  ■  t  *.%  „  „  fn\  -  c/r081  .  where  C  and  f 

We  will  assume  that-  P-^C  ~  pa*F'  '  | 

,  •d/'V-N  ^  p,f^  for  the  law  of  ■ 

|  (jt.  are  certain  constants,  «u-u  ..v\  •}  _  . 


f  „  ovn  ob f p  in  the  p Xs e B s u x* e  0*1^ 

I  motion  of  the  from*  Tncn  ouuxb  ,  * 

|  trlbution  along  Ok  in  the  fonat 

•  i. 


n, 

*  +  r*S£ 


P  (0,  2« 


(2,2) 


j  From  the  conservation  equations  it  follows  that  f  *  */  j 

\/{2*&~)  on  the  front  at  the  point  A.  'Host  61  “?■***  lD*er.  j 
j  eat  centers  about  the  behavior  of  the  solution  when  .  ■ 

S  z— *at.  When  this  happens  it  is  clear  tliac  r  0  >  G  ■  j 

!  Evaluating  the  indeterminate  form  in  (2.2)  «  is  ***!  '-°  j 
1  convince  ourselves  that,  for  example.  Kith  oC=  3,  (the  j 
j  pressure  at  z  =  at)  (at  the  front  of  the  sound  wave)  j 

!  becomes  negatively  infinite.  Thus  in  those  oases  where  j 
px(r,t)  has  a  singularity  close  to  the  front  of  the  j 

j  sound  wave  as  t  —  0,  the  linear  theory  ceases  to  apply,  j 
I  Calculations  show  that  even  when  this  singularity  is  an-  j 
j  sent,  negative  pressures  are  possible  below  the  surface.  J 
j  But  this  appears  to  have  a  definite  physical  meaning.  J 

J  calculation  also  shows  that  tones  of  strong  rarefaction  j 
I  occur  when  the  pressure  on  the  free  surface  drops  sharply  j 
|  with  time,  and  that  they  are  distributed  between  the  free  | 
!  surface  of  the  liquid  and  regions  of  higher  pressure 
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close  to  the  sound-wave  front.  This  result  is  somewhat 

unexpected  If  it  is  kept  in  mind  that  when  fthd  _ 

rT "■  are  constant  in'  time,,  the  pressure  is  essentially 

monotonic  and  positive.;  here;  for  example,  when  VQ  =  a 

the  pressure  changes  according  to  the  linear  equation? 


. "ih  "  at  ' 


(2,3) 


l'  In  order  to  illustrate  the  above  considerations,  j 

I  '  '  ;  f 

1  we  will  -  investigate  the  following  example.  Let  the  | 

I  oressure  distribution  on  the  surface  behind  the  front  be  j 
|  *•  ■*  $ 
|  constant  while  the  pressure  on  the  front  changes  according* 

{  .  .  .  '  .  i 

|  to  the  law  illustrated  in  Fig,  3-  ; 

p,’h\ 


Ha-r  -i 


i 


Ig.  -3 


|  Here  "f  is  some  instant  of  time.  Let  the  pressure  pc 
|  correspond  to  a  front  velocity  Vq  *  a  and  the  pressure 
j  Pg  to  Vg  <  a.  .  ■  •  _ _ 


11 


M,1  J 

In  accordance  with  {2.1}  and  the  law  adopted  for  p. (r,fc) ,  | 

J.,  > 

there  results j 


',  2,  « 


?rnw 


<2,4) 


1^S*+a 


Clearly  for 't  <  t  we  have  t ’  that  much  smaller  than  t 
and  all'  integration  intervals  in  (2.4)  lie  to  the  left 


of  t!  ~  t  in  Flcr 


such  values  of  t  Eq,  (2,3)  holds’. 


We  will  investigate  two  different  cases  for  t  4  it.  j 

1)  z  >  a(t~t) }  here  tl  <  T  also  and  Eq,  (2.3)  holds ;  ; 

2)  z  <  a( t-t) ;  here  t-z/a  >  t,  but  for  some  interval  on  j 

/-*•  ■  •f 

the  z-axis  we  have  t~  ~\|  vQc  +  z^/a  <  13*.  ,| 

t 

Then  (2.4)  takes  the  forms  I 


p(sdt'~~u  ±  \  pu  A 


V 


(2,5) 


At  the  point  z  ®  a(t~t)  there  is  a  pressure  diseonti 


nuit 


p(0t  z,  t)^  ~.pa  f  p*+p 


and  we  take  t  t  and  p0 «  pn,  then  it  is  clear  that  | 

i 

p(OjZj,t)  =  p.  Thus,  the  step-function  pressure  distribu-  ‘j 
tion  .on  the  free  surface  produces  a  nonmonotonic  pressure  jj 

-  1  l 

distribution  differing  from  (2.3)  below  the  surface.,  and  j 


"A*  ?ti*«j**-**U**'  -U  4  t'Atp  **  *  ..  jm*.  4r**lft£4f  *** 


*#*«**•. 


1 


in  certain  cases  produces  zones  of  rarefaction. 
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